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We study the screening of the Coulomb interaction in a quasi one-dimensional superconductor 
given by the presence of either a one- or a two-dimensional non-interacting electron gas. To that end, 
we derive an effective low-energy phase-only action, which amounts to treating the Coulomb and 
superconducting correlations in the random-phase approximation. We concentrate on the study of 
dissipation effects in the superconductor, induced by the effect of Coulomb coupling to the diffusive 
density-modes in the metal, and study its consequences on the static and dynamic conductivity. Our 
results point towards the importance of the dimensionality of the screening metal in the behavior of 
the superconducting plasma mode of the wire at low energies. In absence of topological defects, and 
when the screening is given by a one-dimensional electron gas, the superconducting plasma mode 
is completely damped in the limit g — >■ 0, and consequently superconductivity is lost in the wire. 
In contrast, we recover a Drude-response in the conductivity when the screening is provided by a 
two-dimensional electron gas. 
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I. INTRODUCTION 

The environment has profound effects on the prop- 
erties of quantum systems^. In the case of supercon- 
ductors, it was predicted more than 25 years ago that 
a resistively shunted Josephson junction would expe- 
rience a superconductor-normal transition as a func- 
tion of Rs/Rq, where Rs is the shunt resistance of 
the junction and Rq = h/Ae^ w 6.45 k^l is the 
quantum of resistance^ More recently, a variety 
of superconducting systems, including granulaJ^ , or 
homogeneous^ films, 2D Josephson junctions array^, 
out-of-equilibrium Josephson junction^ and high tem- 
perature superconductor^ were shown to undergo a 
superconductor-insulator transition as the characteristic 
resistance of the system in the normal state increases 
through a critical value on the order of Rq. In those 
cases, the dissipative environment corresponds to the 
measurement circuits or the intrinsic component of nor- 
mal electrons in the system. 

In contrast, isolated superconducting wires with lat- 
eral dimension tq ^ where is the bulk coherence 
length, do not present significant dissipation sources at 
low temperatures. The low-energy modes in an ideally 
isolated superconducting wire are the one-dimensional 
propagating plasmon modes along the axid^. Contrary 
to bulk superconductors, where the plasmon has an en- 
ergy uip^ — \J A-nrLse^ /ni (where Ug is the superfluid den- 
sity and m is the electron mass), in the restricted ID 
geometry of the wire, the long-range Coulomb interac- 
tion is not completely screened and consequently charge 
fluctuations are not shifted to flnite energies in the limit 
g — > 0. The result is a sound-like dispersion relation 
ijp [q) ^ In (l/qro), where the logarithmic factor is a 
remnant of the long-range Coulomb interactions. 

Because of the gapless dispersion relation, quantum 
fluctuations are expected to show critical behavioi'ii'. 



a feature that has attracted the attention of several 
theoretic aP^^^^ and experimentaPsHUl research groups. 




Figure 1: Representation of the capacitively coupled super- 
conducting wire-normal metal system. The metal placed at a 
distance d screens the long-range Coulomb interaction in the 
superconducting wire. In (a) the metal is a diffusive wire, and 
in (b) we consider a diffusive 2D electron gas . 

How this picture (i.e., sound-like dispersion relation 
and critical behavior) is modifled when the coupling 
to the electromagnetic environment is taken into ac- 
count? Intuitively, the presence of a metal at a dis- 
tance d should screen the long-range Coulomb interaction 
for density fluctuations with wavelength q <C , re- 
sulting in enhanced superconducting correlation^^. On 
the other hand, in capacitively-coupled superconductor- 
normal systems, the presence of dissipation in the normal 
metal is known to produce dissipative order-parameter 
fluctuations^^ ■^^ and, from this point of view, screening 
might also carry detrimental effects to superconductivity. 
Moreover, recent theoretical works on related Luttinger- 
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where 



liquid systems coupled electrostatically to metals predict the complete system 
charge-density wave ( CDW) instabilities caused by the 
dissipative environmentP^EIl 

Therefore, a better understanding of the screen- 
ing effects occurring in superconducting wires and the 
consequences to their superconducting properties is 
needed. This is sue is particularl y rele vant to recent 
theoreticaP^I^lIlI] and experiment a P ^ * ^ ^ -1 works showing 
evidence of stabilization of superconductivity in low di- 
mensional systems due to the presence of tunnelling con- 
tacts with normal metallic leads, which suppress of fluc- 
tuations of the superconducting order-parameter. 

In this article we study the effects of the screening of 
the Coulomb interaction in a quasi- ID superconductor 
by the presence of a metal nearby (cf. Fig. [T|). To that 
end, we derive a phase-only action of the coupled sys- 
tem valid at low energies, which amounts to performing 
an RPA-approximation of the interacting problem. We 
specify two experimentally relevant geometries, namely: 
a) a ID and b) a 2D electron gas (IDEG and 2DEG, 
respectively) in the diffusive limit. Our results show a 
rich behavior of the ID plasma mode in the wire due to 
screening effects, and point towards the importance of 
the dimensionality of the screening metal. In particular 
in the case of screening provided by a IDEG important 
frictional effects are observed in the superconductor due 
to the capacitive coupling, and in the limit g — )■ and 
T — > phase-coherence is destroyed and the wire shows 
finite resistivity. In contrast, for a wire screened by a 
2DEG, friction and dissipation vanish in the limit g — 0, 
and the wire is well described by the Luttinger liquid 
picture. 

The article is divided as follows: in Sec. |ll] we de- 
rive a general effective phase-only action for the complete 
superconductor-normal system, in Sec. |III| we present an 
analysis of the screening regimes at low energies for both 
the ID and 2D geometries. Sec. |IV| is devoted to the 
study of the dissipative effects in the dynamical conduc- 
tivity a {q,u}) of the wire, and finally in Sec. |v]we sum- 
marize our findings and present a discussion. The details 
of the derivation of the low-energy effective action are 
given in the Appendices |X] and |B] 



II. MODEL 

In this section we derive a general effective model which 
describes a clean superconductor capacitively coupled to 
a diffusive metal. We leave for Sec. |III|th e specific anal- 
ysis of the systems depicted in Fig! [l] representing a 
superconducting wire of length L and lateral dimensions 
''o ^ ?o coupled to a diffusive metal place d at a dista nce 
d. The derivation of the model is standarcP^IMlIIlll and 
here we only sketch the main steps. We refer the reader 
to the Appendix [X] and to the aforementioned references 
for details. 

In the following we use the convention h — kg — 1. 
We begin our description with the microscopic action of 



^ = drJ2 [dr ra.a (dr - ^a.a + dr H, 
Jo „ „J Jo 



(1) 



where /3 — ]p. The Grassmann field i/'a.CT = V'a,CT(r,T) 
describes an electron in the superconductor for a — s 
(normal metal for a = n) with spin projection a at po- 
sition r = {x,y,z) and imaginary-time r. The chemical 
potential /i = /c|/2to is the Fermi energy in the normal 
state, with kp the Fermi wavevector. The Hamiltonian 
H of the systems is 



(2) 



2m 



(3) 



describes a translationally invariant, clean superconduc- 
tor. Since we will not focus on the details of the pair- 
ing mechanism, here we assume a phenomenological lo- 
cal attractive interaction U < which is responsible for 
(s-wave) pairing at T < Tc- 

The normal metal is described by 



2m 



+ V,UFi^„,4' (4) 



where Vi = Vi (r) represents the weak static impurity 
potential which provides a finite resitivity in the metal. 

Finally, the interaction term of the whole system is 
given by 



Hint j dYidY2 Ps (ri) V (ri - r2, 0) ps (r2) + 
+ \ J c'ridr2 pn (ri) v (ri - r2, 0) /5„ (r2) + 
+ / dridr2 ps (ri) v (ri - r2, d) p„ (r2) , (5) 



where we defined the electronic density operators 
/3a (r) = X]<T V'i.o- (r) V'a.o- (r), and where the domain of 
integration of the variables ri and Y2 is constrained to 
the volume of the superconductor (if a = s) and the 
metal (if a — n). The interaction potential v (r, z) is the 
microscopic long-range Coulomb interaction 



v{y,z) 



where is the dielectric constant of the insulating 
medium between the metal and the superconductor. 

The first step in the derivation of an effective low- 
energy model consists in decoupling the interaction terms 



3 



appearing in and Hint by the means of suitable 
Hubbard-Stratonovich transformations (HSTs). The re- 
pulsive Coulomb interaction H-mt is more conveniently 
decoupled by expressing it in terms of the symmetric and 
antisymmetric density operators 



P± (r) = Ps (r) ± Pn (r) • 



(6) 



With this definition, the interaction term [cf. Eq. ([5|] 
compactly writes 

^int ^7;^ dYidY2 py {Yi)Vy {vi - V2) fly (vi) , (7) 



where we have defined 



Vy (r) 



_ w(r,0) + {v)v{x,d) 



(with V ^±). (8) 



The HSTs to decouple the long-range Coulomb and the 
Hubbard [/ < interactions are implemented by intro- 
ducing the HS fields p^, (r, r) in the particle-hole channel, 
and A (r, r) , A (r, r) in the particle-particle channel, re- 
spectively (cf. Appendix [A|) . 

The next step in our derivation is to introduce an ex- 
tra HS field py (r, r) in order to decouple the quadratic 
term in p,^{r,T), appearing in Eq. (Al). Then, it is 
easy to show that the field p^, (r, t) can be formally 
integrated out, yielding a functional-delta function^*^ 
S [pn (r, r) — py (r, r)]. As noted by De Palo et alW^, this 
fact allows to interpret the new HS fields pi, (r, r) as the 
physical density of the problem, expressed in our case 
in terms of the symmetric and antisymmetric collective 
modes. 

At sufficiently low energies, amplitude fiuctuations of 
the order parameter A (r, t) can be neglected, allowing 
to write A(r,T) — Aoe'^^'"''^-' , with a real constant Aq. 
The phase field 6 (r, r) can be absorbed by a unitary 
transformation of the fermionic field 

(r, r) ^ r) = (r, r) e^^(^^^^'\ 

The derivation of the effective model proceeds with 
the integration of the fermionic fields V'a.o-i ^^id by ex- 
panding the resulting bosonic action around the saddle- 
point in terms of the derivatives of 9 (r, r) and the den- 
sity fiuctuations 5pi, (r, r) , (5p^ (r,T) [cf. Eqs. (AlOl and 
( All l[. This expansion amounts to perfo rniing an RPA- 
approximation of the interacting problerrPSmT 

The last step is to integrate the auxiliary field 
dpi, (r, r), which in the original representation of the den- 
sity in terms Sp^ (r, r) , 5pn (r, t) yields 



ScS = j dvdr -drO (r, t) ps (r, r) + 
1 f ^ 

+ 2/11 (''I ~ ''2, Tl - T2) X 

i=l 

(ri,Ti)y (ri -r2,Ti -r2)<5p(r2,T2) , (9) 



where T) (r, r) is the phase stiffness of the superconductor 
[cf. Eq. ( |A15t [ and 



[Xo.s (r,T)]~' 




5{t) 



w(r,0) v{r,d) 
V (r, d) V (r, 0) 



where 



have 

„2k.r— i 



used the notation [xo,a (r, ''')] 
'''"''Xo i (l^''^ni)' where xo,a{K^m) 



the bare density-density correlator (i. e., ob taine d with 



the Hamiltonian defined in Eqs. ( |A12[ ) and ( |A13 | 
Here we have used the notation in Fourier space (k,w,„) 
with k the momentum and uj„i = ^jf^ the bosonic Mat- 
subara frequencies^. 

Note that at T = and in absence of quasiparticle exci- 
tations, the whole electronic density in the superconduc- 
tor corresponds to the superfiuid density. Consequently, 
the field 5ps (r, r) physically represents the fiuctuation of 
the Cooper-pair density at point (r, t). 

An interesting aspect of the effective action in 
Eq. ([9]) is that the first term (i.e., coupling be- 
tween the total density of Cooper-pairs ps (r, r) and the 
phase-field 0(r, r)) appears naturally as a consequence 
of the well-known number-phase commutation-relation 
[ps {r) ,6 {r')] = iS (r — r') occurring in the supercon- 
ducting groundstate^. 

Besides the contribution of soft modes, encoded in Eq. 
([9]), in low-dimensional superconductors there are also 
stable topological excitations which contribute to the ef- 
fective action. These are the well-known classical vortex 
(in 2D) and the phase shps (in ID) excitationd^^HUl, Fo- 
cusing in the ID case, a phase-slip is a region of size ^0 
where the order parameter temporarily vanishes, allow- 
ing the field 9 (r, r) to perform a jump of ±27m (with 
ri integer) across it. For wires in the limit of very low 
superconducting stiffness, phase slips are an important 
source of momentum-unbinding, and a releva nt contri- 
bution to the action in the RG-sens^^SHIlIl Indeed, 
it is believed that the eventual destruction of the su- 
perconducting state in isolated ultrathin wires occurs 
through the pr oliferation of qu antum phase slips/anti 
phase slips pa,ii^MMMIMM31i^ in what constitutes the 
quantum analog in 1+1 dimensions to the classical 
Berezinskii-Kosterlitz-Thouless (BKT) transition in two 
space dimensions'^^. 

Note that our derivation does not account for the pres- 
ence of phase-slips. Consequently, our results will only 
apply far from the BKT transition and far from the (non- 
superconducting) phase where the effect of phase slips 
dominates. 

In the following we analyze the generic action of Eq. ^ 
for the different configurations of Fig. [T] 
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III. SCREENING REGIMES 

A. Unscreened isolated wire 

Let us first explore the instructive case of a super- 
conducting wire ideally isolated from the environment. 
This situation corresponds to the normal metal placed in- 
finitely far from the superconductor (i.e., d — >■ oo), which 
results in the decoupling of their dynamics. For a very 
narrow superconducting wire with tq <SC ^oj the depen- 
dence of the fields 9 (r, r) , Sps (r, t) on transverse dimen- 
sions can be neglected, reducing to {9 (r, r) , Sps (r, r)} ^• 
{9 (x) , Sps (x)} where the compact notation x == {x, t) 
has been used. The effective action of the superconduct- 
ing wire writes most conveniently in Fourier space 



^ 2^ E '^"^ (q) Ps (q) + 1''^ (q) (q) 



1 + Xo,s (q)w(g,0) 2 

TTT \Ps[<V\ 

xo,s q 



(10) 



where we have used the notation q — {q, —OJm) (with q 
the momentum along the wire), and the property of real 
fields 61* (q) = 6'(-q), (q) = p^-q). The Fourier 
transforms V (q) and xo,s (q) are defined as 

V{q)^ f dT [ dx [ dydze-''^'^V{Y,T), (11) 
Jo Jo J 

Xo,s (q) = f dT ( dx ( dydz e"'''^ xo,s (r, r) . (12) 
Jo Jo J 

At this point is relevant to calculate the Fourier trans- 
form of the Coulomb potential, which can be approxi- 



mated as V (r, z) 



(i.e., only dependent on 



the spatial coordinate x), and cut off at short distances 
by the radius tq. Therefore we have 



v{q,d) = Ko { \q\ 



r^ + d^ 



(13) 



where Kq (C) is the zeroth-order Bessel function, which 
verifies the limit lim^^^o Kq (C) — In ( | j — 7, with 7 



the Euler gamma constanpSl From the above Eq. ( 10 ) 
we obtain the phase-only action in the limit q — >■ 0, by 
integration of the field ps (q) 



n 



2/3L ^ 



In 



(k^ol) 



|e(q)r, (14) 



where Vq = limq_^o (q) = ^ (cf- Appendix 

The minimization of the effective action Eq. ([141) al- 



lows to obtain the equation of motion for the phase- 
field and to recover the dispersion-relation predicted for 



the ID-plasma modd^^l^iJ upon analytical continuation to 
real frequencies iujm — > w + iO~^ 



to' (q) 



0. 



(15) 



Let us now concentrate on the superconducting proper- 
ties of the wire. It is well-known that long-range order of 
the order parameter in ID quantum systems is not possi- 
ble, due to presence of strong quantum fiuctuations and, 
strictly speaking, only quasi-long-range order, character- 
ized by a slowly decreasing order-parameter correlation 
function 



i^(x)^(A* (x)A(O)) 



(16) 



can exisiF^^^. In the case of the isolated wire, the phase- 
correlation function calculated with the effective phase- 
only action Eq. (14) write^^^^ 



Tr[9{^)-9m 



where K 



Voir 



1 



and u 



In 



y/ x'^ + u^t'^ Inr 
ro 



3/2 



(17) 



As compared with 



the case of a ID superconductor with short-ran ge r epul- 
sive interactionJii', the phase correlator of Eq. (17) pro- 



duces a relatively fast decrease of the order-parameter 
correlation function Eq. (16), as a consequence of the 



long-range Coulomb interaction, which is not completely 
screened in the ID geometry. Consequently, density fluc- 
tuations are suppressed in the limit q — )■ (JSl and su- 
perconductivity, which benefits from fiuctuations in the 
density, is suppressed. 

A natural step to take in order to diminish the detri- 
mental effects of the Coulomb interaction in the ID ge- 
ometry, is to screen it by the means of a metal placed 
nearby. This is the subject of the subsequent sections. 



B. Screening by a diffusive metallic wire 

We now concentrate on the system depicted in Fig. 
[ija) . For simplicity, we consider the case of two geomet- 
rically identical cylindrical wires. Extensions to other ID 
geometries are straightforward. We assume that the nor- 
mal metal is only one-dimensional with respect to density 
fluctuations p„ (q) with spatial wavevector q satisfying 
the condition qr^ ^ 1. Note that this condition does 
not necessarily imply that the normal wire is electron- 
ically ID (i.e., it does not imply the existence of only 
one electronic conduction channel). Indeed, in what fol- 
lows we assume a normal metal with a large number of 
channels iVch ~ (feFJ^o)^ ^ 1- This fact, together with 
the additional assumption of a very weak disorder poten- 
tial, allows to neglect Anderson-localization effects (i.e.. 
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Figure 2: Screening regimes for a superconducting wire 
screened by a diffusive fDEG. Tlie curve to — Dq^ separates 
tlie regime of static screening uj <^ Dq^ (light gray area) from 
that of dynamical screening uj 3> Dq^ (white area). The dis- 
persion relation of the ID plasma mode (thick solid line) is 
obtained from the solution of Eq. \2Q\ . For typical experi- 
mental values (cf. Table|l|, the dispersion relation crosses over 
from the static regime to the dynamical regime, and eventu- 
ally the mode is completely damped. The unscreened regime 
of frequencies lo > Dg^^A/"" ^ 
dark gray area. 



In — corresponds to the 



L <C ^wirc, where fwire is the localization length in the 
diffusive normal wire). 

In that case. Eq. (pi reduces to 



""jbiY. (q) p*s (q) + (q) 1^ (q)l' 

q 

-t-pt (q)y(q)p(-q)], 



(18) 



where the subindex q in ST indicates the effective di- 
mensionality of the metal. The integration of the density 
modes ps (q) and pn (q) in the above expression allows 
to obtain the result 



(1) 



r, 1 A 



2(3L 



+ (xo.s (q) + xo.n (q)) V (q, 0) 
xo,s (q)(i + xo,n(q)^'(g,o)) 



Xo.s (q) xo,n (q) (v {q, of - V (q, dy 

Xa,s (q)(l + Xo,n (q)w(g,0)) 



+ q'V{q)U9iq)\ 



(19) 



In the following we focus on the experimentally rele- 
vant regime d « rg ^ q^^. In that case the quantity 

V {q, of - V (g, df - [in ^fc^O drops from Eq. ([l9 ) 
and the expression simplifies to 



QW 

^(1 



(1) 



— T 

q 

..2. 



Xo,5 (q) [1 + Xo,n (q) V (g, 0)] 



4 1 + [Xo,. (q) + Xo.n (q)] v {q, 0) 



+q'V{ci)}\0{q)\' 



(20) 
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Figure 3: Real and imaginary components of the ID plasma 
mode cj (g) , obtained from the equation of motion of the ac- 
tion Eq. ( |20| . The real part (solid line) gives the dispersion 
relation, while the imaginary part (dotted line) represents the 
damping of the mode. As in Fig. [2] the curves have been cal- 
culated for realistic experimental parameters (cf. Table |l| . 
The curve Dq^ (blue dashed line) is shown as a reference. 



For a weakly-disordered diffusive electron gas with 
elastic mean-free path l^. and scattering time Tg = le/vp, 
where is the Fermi velocity, the disorder-averaged 
density-density correlation function [cf. Eq. (A13jl at 



energies \ujm\ < f^nd momentum q < writes 



XO,n 



(q) ^ ^K.i 



Dq^ 



D 



Dq^ 



(21) 



where JV^^ uj is the ID density of states at the Fermi 
energy in the normal metal, and D = 'e/''"e is the diffu- 
sion constant in ID. The factor 2 accounts for the spin 
degeneracy. 

Note that the susceptibility Xo,™ (q) [cf. Eq. (21l[ is 
non-analytical in the limit q — > for a normal diffusive 
metal. On the contrary, for the superconductor the pres- 
ence of a gap in the excitation spectrum allows to obtain a 
well-defined limit hmq_j.o Xo,s (q) — 7A/'° i£), where A/'f i^, 
is the linear density of states in the superconductor (in 
the normal state) at the Fermi level, and 7 is a numerical 
coefficient of order 1 [cf. Eq. (Bl)[. 



The plasma mode obtained from the equations of mo- 
tion derived from Eq. ( 20 1 is plotted in Fig. [2] (thick 
solid line). Due to the complexity of the screening pro- 
vided by the diffusive IDEG, it is instructive to derive 
analytical expressions valid in the limiting cases of static 



(i.e., \u!m\ ^ Dq^) and dynamical (i.e. 
screening. 



Dq^ < \0Jrn\) 



1. Static screening limit \(jJm\ ^ Dq'^ 



This limit corresponds to the region \uj„i\ ^ Dq^ (see 
light gray area in Fig. [2]). In this case, the susceptibility 
in the normal metal can be approximated as xo.n (q) 

2AC? (1 ~ ~ 



[cf. Eq. (21 1[ . Then Eq. (20) can be 
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written as 
1 



2/3i ^ 
q 



Xs (0)^^ 



a. |a;„ 



Dg2 



|e(q)r 

(22) 



with Xs (0) 



susceptibility of the wire and a 
that in the Hniit a 



the effective static RPA- 



Note 



0, the above action corresponds to 
a Luttinger Hquid action with short-range interactions^^. 
In the case of a ID geometry of Fig. [l|a), the screening 
length is given by the distance d. 

In the more general case of a > 0, the term ^ "j^"^ 
introduces dissipation in the plasmon mode. From Eq. 



( 22 1 , the dispersion relation for the plasma- mode writes 
au (q) ' 



Xs (0)' 



= 0. 



(23) 



This equation holds provided the consistency condition 
|(u;(q)| ^ Dq^ is verified (cf. solid line in Fig. [2|. In 
Fig. [3] we show the solution of the above Eq. ( |23| as 
a function of q. Note that while Re[cj (<?)] follows an 
approximately linear dispersion relation, the imaginary 
part takes a constant value in the regime |w ((?)| ^ Dq^, 
meaning that the plasmon mode acquires a finite width, 
which in the perturbative limit a — >■ writes F (q) = 
-Ini[uiq)]^^^ (cf. Fig.|3). 



represents the effect of quantum fluctuations induced by 
the charging energy of the superconducting islancPSl. As a 
result, dissipation turns out to be beneficial to supercon- 
ductivity, through the quenching of phase fluctuation^^. 
However, in our case, the form of the action in Eq. 



(25 1 is qualitatively different, since the term ^ w„j is ab- 
sent from the action (actually, it is the dynamical term 
itself which becomes a contribution ~ q^ |wm|). This has 
detrimental consequences for the superconductivity in 
the wire, as can be seen directly from the equation of mo- 
tion for the fleld 6, which gives uj (q) ~ —iDo/ {DAf^ , 
indicating that the original plasma mode is completely 
damped and vanishes in the limit q ^ (see Fig. |2|. 
Indeed, expressing the action Eq. ( 18 ) in terms of the 
dual fleldSI ^ (x) , deflned as 



Sps (X) EE V0 (x) , 



we obtain the equivalent description 



(26) 



OW 

-^(1 



(1) 



1 1 

2^/3L 




+9' 



1 + [Xo,s (q) + Xo,» (q)] v {q, 0) 

xo,s (q) [i + xo.n(q)f (g,0)] 



In the regime of Eq. ( 24 1 , we can approximate the action 

by 



2. Dynamic screening limit \ujm\ 3> Dq^ 

For realisitic estimates of the experimental parame- 
ters (cf. Table |l|, our results indicate that the regime 
\ujm\ ^ Dq^ (white area in Fig. [2]) is the most relevant 
for experimental studies on today's accesible wirei^SHUl, 
Replacing Eq. (21 1 into Eq. (20) we note that if the 
condition 



2 

Dq'' « \lu„,\ « Af^^ ^j,Dq^ In , 



(24) 



is fulfilled, the action in Eq. ( 20 1 can be approximated 
as 



The action Eq. (25 1 indicates that phase fluctua- 



tions show dissipative dynamics (encoded in the term 
~ g2 loi^l) as a consequence of the coupling to the dissi- 
pative processes in the IDEG. In other words, the super- 
conductor "inherits" the dissipation in the IDEG through 
the Coulomb interaction. 

Note that a term ^ q^ |a;„i| has been studied in the 
context o f resistiv ely shunted Josephson junctions arrays 
(RSJJAs)ISE2llll. In that case, the term ~ q^ ap- 
pears in addition to the dynamical term ^ w^, which 



S\ 



1 1 



~ 27r2 PL 
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7AA% + 2AA„%|3il„(,,o) 

^ N0(q)l , 



,2 



2^2 /3L { 2U{['^D + 4Po ^ iKu '"^ ^"^^'^ 



ID 



q^ Icj,, 



in Eq. ( 25 I trans 



which shows that the term 
lates into a relevant term ~ jw^l (in the RG-sense) 
when expressed in terms of 4> (q) ■ Another way to see 
this detrimental effect is through the order-parameter 



1 cor 
^a) ^ ^ E WliDDq^ M + ^09^] \e (q)l' . (25) Eq 



correlation function F (x) = AqC 



T,[e(x)-e(o)]=^ 



[cf. 



(I61], which vanishes due to the infrared di- 

2\ 



vergence of the phase-correlator (T,- [6* (x) — [Q)]' 

1 — cos q.x 



— ?► 00. 



(cf. 



PL 2^I°Dq^\l^J„^\+Vaq^ 

Only at high-frequencies ^ 
dark gray area in Fig. |2|, and provided Eq, 
still valid, or in the limit of very low electronic density of 
states in the IDEG, we recover the action of Eq. ( 14 ) de- 



(211 is 



scribing again unscreened plasma modes. Physically, at 
such high frequencies the response xo,n (q) of the IDEG 
vanishes and the superconducting wire is effectively un- 
screened. 
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ro ~ d 


L 




A/-0 ~ A/-0 

-'\'s,lD — ■'^n,lD 


Ao 


D 






1.20 




10 nm 


100 /im 


8-6 10^^ ^ 


10^9 1 

m. J 


1 K 


0.01 


1 


100 nm 


1 nm~^ 


10^'^ 



Table I: Parameters used in the calculations. Order-of-magnitude estimations of ro and L have been extracted from experiments 
on superconducting aluminium wires with coherence length estimated as ~ 100 nm [cf. Ref. I18| . 



C. Screening by a diffusive metallic film 

Now we focus our attention on the system of Fig. [ijb), 
which represents a superconducting wire coupled to a 
normal diffusive film of width Wfum- In this case, the 
effective action Eq. (pi writes in Fourier space 



2«iz7 E| [^O.n (q, k) + V (q, k, 0)] |p„ (q, fc)|' + 



+ V (g, fc, d) [pi (q) p„ (q, k) + ps (q) p,*, (q, k)] 



(27) 



where the Coulomb interaction is [compare to Eq. ( 13 1] 

v(q,k,d) = , „ ,„ ■ 

In this case, the presence of the superconducting wire 
breaks the translational symmetry of the system in the 
direction perpendicular to the wire. Consequently, the 
perpendicular momentum k in the plane is not conserved 
and the Coulomb interaction couples the density modes 
in the wire with momentum q to all the modes in the 
plane with momentum k. As before, in order to obtain 
an effective model for the phase field, we must integrate 
over the density fields p„ (q, k) and ps (q), which yields 



(2) 



1 1 

2/3Z 



E- 



xo.s (q) 



4 1 + xo,s (q) [v (g, 0) - v^s (q)] 



+ 9'2?(q) |e(q)r 



(28) 



where VeS (q) is an effective ID-potential encoding all the 
screening properties of the diffusive film 



Vcs (q) 



1 V (g, k, d) xo,n (q, fc) 

L± ^ i + w(g,fc,o)xo,n (q,fc)' 



r 



In this case, the susceptibility at low energies write^^ 

D (g2 + p) 



xo,n (q, k) 



o/v/'O 



(30) 



'Diq^ + k^) + \LU„,\' 

where Af^ 2D i^ the 2D density of states at the Fermi 
energy in the normal metal. Here again, we neglect 
Anderson-localization effects in the metal by assuming 
that the length of the wire is L <C &im, where ^mm is the 
localization length in the film. 

In what follows, we derive analytical expressions for 
the model in the limiting cases of static and dynamic 
regimes. 



1. Static screening limit |aj„j| Dq^ 

This region corresponds to the ligh gray area in Fig. 
([4]). A relevant length scale which naturally appears is 
the 2D Thomas-Fermi screening length, A^p = 2e^Af" — ' 
beyond which the long-range Coulomb potential is com- 
pletely screened. This quantity defines the Thomas- 



Fermi wavevector /ctf = 
tally relevant limit kq^pd 3> 1, 
Woff (q) reduces to 



In the experimen- 
the effective potential 



■^cff (q) 



2e2 



A'o {2qd) - 



2 Dqh 



TF 



(31) 



The first and second term in the above expression are 
consistent with the static and dissipative contributions, 
respectively, to the effective screened interaction ob- 
tained for a Tomonaga-Luttinger liquid electrostatically 
coupled to a diffusive 2DEG (cf. Ref. The 
static screening provided by the 2DEG [first term in Eq. 
(31 1[ cuts the logarithmic divergence of the bare intrawire 
Coulomb interaction v{q,0). The relationship between 
the second term in Eq. ( |31[ ) and the dissipative contri- 
bution ~ q \ujm\ in Ref. 23 can be made explicit with the 
introduction of the field (j> (x) , defined in Eq. ( 26 1 . 

In the limit qd ^ (with \uJm\ ^ Dq^), the effective 
potential Eq. (31) can be further simplified to Veg{q_ 



2e 



2 DqkTF 



and when replaced in Eq. ( 28 1 



yields 



(2) 



1 

2j3L 



q 



'-Xs (0) 



9 9 

T^e/Xs 



(0) 



Dk^FQ 



(29) where Xs (0) Xo,s (0) [l + ^Xo^s (0) In (2d/ro) 
the RPA susceptibility of the wire. 
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Note that in the limit £> — >■ cxo (no dissipation in 
the normal metallic film), we recover again the action 
of a Tomonaga-Luttinger liquid with short-range inter- 
actions j_with plasma modes obeying a linear dipersion 
relation^. In the more general case of a finite D, we 
have the relation dispersion 



(q) Xs (0) + i 



TTUJ^ {q)xl (0) 



Dqh 



TF 



0, 



(32) 



which describes plasma modes with approximately lin- 
ear dispersion relation, and with a width T {q) ~ 



-Im \uj (gr)] 



1^0 

2e,DfcTFXs.o(0) 



q. Note that this result only 



applies in the limit |w {q) \ <^ ^Dq^, Ii^YTSt}' ™^ even- 
tually breaks down in the limit q 0, meaning that this 
is not the relevant regime at low energies. 

2. Dynamical screening limit \LUm\ ^ Dq^ 

With realistic estimates for the experimental parame- 
ters (cf. Table [l]), the regime jw^l ^ Dq^ is the most 
relevant in practical realizations. In this regime (white 
area in Fig. |4]) , the effective potential v^s (q) in Eq. ( 29 ) 
can be approximated as 



VeB (q) 



2e2 



/ (2fcTFd) - / 



2 \ujm\ d 



Dh 



TP 



where we have defined f {z) = — e^Ei(— z), with Ei (a;) 
the exponential integral function '^^ . If the additional con- 
dition \(jJm\ ^ ^^d^ holds, the effective potential can be 
further simplified to 

2e2 /2KrJ , 
i;efr(q) ~ In -jrr—d 



Using this expression and Eq. ( 13 1, the phase-only action 
of Eq. (l28l) writes 




Xo,s (0) 



2e2 



q^vA\e{ci)\\ 



resulting in the equation of motion (in the limit q 
" ,2e2 



UJ 











droquj 





= 0. 



^0) 
(33) 



In Fig. |4] we show (solid line) the dispersion relation 
obtained from Eq. (33l (i.e., real component of uj{q)). 
Contrarily to the case studied in Sec. IIIB the resulting 
plasma mode is not damped in the limit g — >■ 0, i.e., a dis- 
persive real component survives. In order to investigate 
the dynamics of the phase-field at low-energies, we study 




0.0001 0.0002 0.0003 0.0004 



q [nm '] 



Figure 4: Screening regimes for a superconducting wire ca- 
pacitively coupled to a diffusive 2D EG . The dispersion rela- 
tion (solid line) results from Eq. (|33|l, valid in the regime 



Dq^ < |a;„| «; As in Fig. [2] the dashed line uj = Dq 

separates the regime of static (light gray area) from that of dy- 
namic (white area) screening. For the same parameters as in 
Fig. |2} the pi asma mode crosses over between the two regimes 
at higher frequencies. The unscreened region corresponds to 
the dark gray area. 



the ratio Im (q)] /Re [cj (g)] vs. q (cf. Fig. [s]). From 
Eq. ( 33 I it is possible to show that in the limit q ^ Q 



Im [uj {q)] TT 
Re [cj {q)] ~ 4 



In 



Dh 



TP 



^q'dro^ 



Vn 



(34) 



meaning that the width of the plasma mode decreases at 
low energies, resulting in a well-defined excitation. Note 
the difference with respect to the screening provided by 
a IDEG, where the damping of the plasmon was com- 
plete in the limit q — > 0. The origin of this difference 
lies in the additional degree of freedom k (momentum 
perpendicular to the wire), which smears (upon integra- 
tion) the dependence on the damping factor |w,„| in the 
susceptibility [cf. Eq. (30)] . 



lmMq)yReMq)] 
k 

0.05 
0.04 
0.03 
0.02 
0.01 



0.0001 0.0002 0.0003 0.0004 



q [nm '] 



Figure 5: Ratio vs q for typical experimental param- 
eters (cf. Table The plasma mode becomes better defined 
as g — >■ 0. 

In the regime of frequencies ^'^^^^ <C <^ Dk^p, 
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the effective potential Wcff (q) writes 



VcB (q) 



2e' 



1 



and the phase only action is 



LOr, 



(2) 



1 1 



Xo,. (0) 



+ ¥xo.(0) [in(^) 



In this limit, the equation of motion of the field 6 (q) is 



o2 r 



In 



iDh 



TF 



qrn J 2duj (q) 



Note that in this regime, the dissipative effects are even 
weaker and the dispersion relation ressembles that of the 
(unscreened) Mooij-Schon mode. Eventually in the limit 
^ Dk^p, the response xo,n (q) of the 2DEG van- 
ishes and the wire is effectively in the unscreened regime 
where the Mooij-Schon plasma mode of Eq. ( 15 1 is fully 
recovered. 



IV. DISSIPATIVE EFFECTS IN THE DYNAMIC 
CONDUCTIVITY 

In this section we study the consequences of thedissi- 
pative effects on the dynamic conductivity of the wire 
a{q,uj), i.e., the ratio between the current density and 
the local electric field j{q,uj) — a {q,uj) E {q,uj). This 
quantity is of interest because its real part Re [a {q, cj)] 
provides information on the dissipation and absorption 
properties, which result in our case from the friction me- 
diated by the Coulomb interaction'^. 

The response of the system to an external electromag- 
netic field can be obtained by the means of the minimal 
coupling — iV — > — iV — in the microscopic Hamilto- 
nian Eq. ([s]). For a superconducting wire at T = and 
in absence of quasiparticle excitations, the total current 
density is given by 

j (x) = jp (x) + jd (x) , 
jp (x) = -^VoWe (x) , 

Jd(x) = -(^^^ 2?oA(x), 

where jp and jd are, respectively, the paramagnetic and 
diamagnetic contributions to the current density. The 
linear response to an applied electromagnetic field is 
given by the current-current susceptibility of the wire 



X3J (q) 



SlnZ 



Op(q)ip(-q)> -2^0 



A=0 




Figure 6: Dynamic conductivity Re a (g, uu) of a supercon- 
ducting wire dynamically screened by a diffusive IDEG. The 
plasma mode, which is better defined at high energy and mo- 
mentum becomes completely damped in the limit {oj, g} — >■ 
by the effects of the dissipative environment. 



The conductivity is in turn related to th e current- 
current susceptibility by the relatiorfSEIl a (q) — 

— ^'^^"^"^ , upon analytical continuation to real frequencies 

i'i)\iLo -^u+iS = '^(li^)- III terms of the phase field 
0{q), the conductivity reads 



a(q)^- 



2 r 



73)2 „2 

— (^(q)^(-q)) 



(35) 



Let us first study the response of an ideally isolated 



wire (cf. Sec. Ill A I to the electromagnetic field. At 
T = we obtain 



Re [a {q, uj)] 



/2e 



(36) 



where uj (q) = \ ^ In ( j^f^) Q^T^o is the energy of the 



Mooij-Schon plasmon [cf. Eq. (15)]. The real part of 
the conductivity tells us that the system absorbs en- 
ergy at the frequency w = '^(9)1 which in this case 
are well-defined excitations (i.e., delta- functions). Note 
that in the limit q 0, Eq. (36 1 allows to recover 
the Drude pe ak at a; = 0, which is expected for a 
superconductoi^JJ^. 

Let us now consider the case of a wire in the proximity 
to an electron gas. We first study the case of screening 
by a diffusive IDEG (cf. Sec. |III B2[ ), where the effects 
of dissipation are at their strongest. Using the action of 



Eq. ( 20 1 to evaluate the formula of the conductivity [cf . 
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Re [criq=0,aj)/o-dc] 




(jj 



8 10 oj 







Figure 7: dc-conductivity a (uj) = Re [a {q = 0,aj)] of a super- 
conducting wire screened by a diffusive IDEG. Tfie values on 
tlie axis are normlized to a^c ~ cr(tj = 0) — (^)^ 2DAfn,iD 
and ujo = 2j\f°'° d ' plasma mode at finite q, is com- 

pletely damped in the limit q — (cf. Fig. [6]|. 



Eq. (35l] we obtain the expression (valid at T = 0) 



2 „2 



■Ini 



2_ (uj + iO+) 



XO^S (0) [l+xan(g,^)«(g,0)] 
l+[xo,s{0)+xh% {q,oj)]v{q,0) 



sequences in the excitation spectrum of the ID supercon- 
ductor. 

As we mentioned before, far from the BKT quantum 
critical point, phase slips are an irrelevant perturbation 
(in the RG-sense) . In the case of Luttinger liquids with 
short-range interactions, the perturbative effect of phase- 
slips generates a power-law resistivity g ~ T'^, with z/ 
a positive exponent Although we have neglected the 
perturbative effect of topological excitations in our for- 
malism, the fact that a finite resistivity at T = appears 
in the superconducting wire indicates that their effect 
in the conductivity might be negligible as compared to 
those induced by dissipation in the electron gas. 

Re o-(q,aj) 



(37) 




where xh%il^^) = l™5^o+ [Xo,« (q)],^,,^^^+,5 is the 
(disorder-averaged) retarded density-density correlation 
function in the IDEG. In Fig. |6]we show the result for 
Re [a (9, w)] of Eq. (37 1 in the plane q — oj. The disper- 
sion relation Re [uj {q)\ vs. q (thick solid line in the bot- 
tom plane) was calculated numerically from Eq. ( 20 I and 
corresponds to the same plot of Fig. |2] As mentioned be- 
fore, the absorption peaks of Re [a {q, u)] are centered at 
the frequency Re [uj (q)] of the plasma mode. The curve 
Dg^ (blue dotted line) is plotted in the bottom q — uj 
plane to visualize the different screening regimes. Note 
that the dissipative effects in the normal wire (encoded in 
a finite value of the diffussion constant D) are manifested 
in this figure through the finite width F (g) ~ — Im [cj (q)] 
of the plasmon peaks. Note in addition that the constant 
width F (q) in the regime \uj {q)\ ^ Dq^ is consistent with 
the result for Im [uj (q)] of Fig. [3] 

As q — > 0, the plasmon peak merges smoothly into 
the dc-conductivity value a^c = (^) "^DM^u^. which 
exactly corresponds to the dc-conductivity of the IDEG 
(cf. Fig. |7]). This expression is obtained by replacing 
the expression of the action of Eq. ( 25 1 in the general 
expression of the conductivity Eq. (351. Physically, this 



means that the Coulomb interaction produces friction in 
the superconductor through the dissipation existing the 
IDEG. It also indicates that the original plasma mode 
is no longer a well-defined excitation of the system, and 
that the electromagnetic environment have profound con- 



Figure 8: Dynamic conductivity Re [cr {q, uj)] of a supercon- 
ducting wire dynamically screened by a diffusive 2DEG. The 
plasma mode is completely damped at high energy and mo- 
mentum, but in the limit {uj, g} — >■ the effects of dissipation 
vanish. 

In the case of screening by a diffusive 2DEG, the ex- 
pression of the conductivity is given by the expression 
(vahd at T = 0) 



Re [a {q,uj)]^Vl 



2eVq^ 



Im 



2^ (uj + iO+) 
q^Vo - ^ 



Xo.s (0) 



1 + Xo,. (0) [v{q,0)~vl^Hq,Uj)] 



where vl"^ {q,uj) = \im.s^o+ [veS {q)],^^^+is- Our main 
results in this case are presented in Fig. [8| Contrarily to 
the case of Fig. [6] the plasmon peaks are worse defined at 
high energies, while at low energies the width of the peak 
centered at Re [uj (g)] decreases and eventually vanishes 
in the limit g — > 0, in agreement with Eq. (34 1 and 
Fig. Eventually, the plasmon peak merges into the 
superconducting Drude peak at w = 0. 

The presence of an additional degree of freedom (i.e., 
momentum in the plane perpendicular to the wire) is of 
central importance to understand the vanishing of dissi- 
pation. Indeed, even in the dynamical screening regime 
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Dq^ <^ \oJm\ for which one would naively think that dis- 
sipation effects are dominant, the presence of perpendic- 
ular wavevectors k satisfying the condition \oJm\ ^ Dk^ 
make the dissipative processes less important. Note in 
addition that this condition is more easily satisfied in the 
limit \ujm\ 0. These qualitative phase-space considera- 
tions allow to understand the behavior of the effective ID 
potential v^s (q) of Eq. ( 29 1 , which produces a weaker 
(i.e., logarithmic) dependence on the term encoding 
the dissipation. The net result is that the ID plasma 
modes are better defined in the limit g — ^ and the fric- 
tional effects vanishes strictly in the thermodynamical 
limit L — > cxD. 



V. DISCUSSION AND SUMMARY 

In this article we have studied the effects of the lo- 
cal electromagnetic environment, provided by the pres- 
ence of a non-interacting electron gas, on the low-energy 
physics of a superconducting wire. In particular, we have 
focused on the derivation of an effective phase-only ac- 
tion, starting from the microscopic Hamiltonian of the 
system. We make extensive use of the path-integral for- 
malism, which enables to decouple the superconduct- 
ing and long-range Coulomb interactions by the means 
of Hubbard-Stratonovich fields, and to expand the re- 
sulting action in terms of quadratic deviations of these 
fields around the saddle-point (i.e., Gaussian fiuctua- 
tions) . This treatment is equivalent to performing the so- 
called RPA-approximation of the interacting problenJ^. 
We have studied two cases in particular, namely, the 
screening provided by a diffusive IDEG, and a diffu- 
sive 2DEG placed at a distance d ~ tq from the wire, 
with To its radius. This would be the relevant situation 
in practical realizations in, e.g., superconductor/normal 
heterostructures made by the means of the ferroelectric 
field-effect in Nb-doped SrTiOg layerj^, or in electrically 
controlled LaAlOa/SrTiOa interfacetPlkl. 

It is of interest to put our results in the context of 
other works dealing with electrostatically coupled ID sys- 
tems. Among these, the Coulomb drag effect'"'", where 
a finite current Ii is driven in one (the "active") sys- 
tem, and a finite voltage V2 is induced in the other 
("passive" system), has received a great deal of atten- 
tion both theoretically'^^ and experimentalljlSSHSU^ Al- 
though closely related, the focus of our work is on the 
equilibrium properties of the wire. 

From the theoretical point of view, our work differs 
from the usual Tomonaga-Luttinger liquid description of 
a purely ID (i.e., one electronic conduction channel) con- 
ductor, where the main mechanism of momentum trans- 
fer is backscattering^'^ ■'''^ '"'^ . Indeed, it is worth to note 
that intra- and/or interwire backscattering effects are ab- 
sent in clean wires with a large number of electronic chan- 
nels, and this fact is correctly reproduced by our effec- 
tive coarse-grained theory [cf. Eq. ([9|]. Therefore, in 
the language of Tomonaga-Luttinger liquid physics, our 



treatment amounts to retaining only forward scattering 
processes. 

Our results point towards a rich behavior of the ID 
plasmon mode in the wire, determined by the diffusive 
modes in the electron gas. Independently of its dimen- 
sionality, in the static screening limit |a;,„| ^ Dq^ , the 
wire has a plasmon excitation which follows approxi- 
mately a linear dispersion relation. One could naively 
think that in that regime dissipative effects are negligi- 
ble. However, the complete solutions of Eqs. (23 1 and 



(32 1 indicate that this is not the case. Indeed, we obtain 
sizable dissipative effects even in the limit jw (g)| <C Dq^, 
which manifests itself in the broadening of the ID plas- 
mon mode (cf. Figs. [6]and[8|. Although technically chal- 
lenging from the experimental point of view, this broad- 
ening could be seen in experiments of resonant inelastic 
Raman light-scattering^ or in optical measurements of 
the dynamic conductivity or the reflection coefficienlF^. 

On the other hand, our results reveal that the dynam- 
ical screening regime Dq^ <^ \'^m\ should be the most 
relevant for experimental realizations (cf. Figs. [2] and 
[4]). This is more or less evident from the fact that the 
plasma mode essentially follows a linear dispersion in the 
limit g — >■ 0, while the boundary between the dynami- 
cal and the static screening regimes (determined by the 
diffusive modes in the electron gas) is ^ Dq^ . More im- 



portantly, in the limit Dq <^ 



the dimensionality 



of the electron gas is of central importance to determine 
the low-energy properties of the wire. If the screening is 
provided by a IDEG, its dissipative processes are more 
efficiently transferred to the superconducting wire in the 
limit q ^ 0. As a consequence, the plasma mode becomes 
an ill-defined excitation and the superconductor shows a 
finite dc-conductivity in the limit a; = (cf. Figs. |6]and 
I?]). This effect could be seen, e.g., in dc-transport exper- 
iments on capacitively coupled superconducting/normal 
wires systems (cf. Fig. [t]). 

When the screening is provided by a 2DEG, acoustic 
plasma modes with a vanishing width are recovered in 
the limit g — >■ 0, which allows to neglect the dissipative 
effects due to the Coulomb interaction with the metal 
(cf. Figs. |5]and[8]). The reason for this lies in the exis- 
tence of the additional degree of freedom in the electron 
gas (perpendicular momentum k) , which produces (upon 
integration) a weakening of dissipation effects. At this 
point it is tempting to speculate that a semi-infinite 3D 
metal, or a superconducting wire embedded in a 3D nor- 
mal matrix, would provide an additional degree of free- 
dom (momentum k' perpendicular to the wire and to k) , 
and would weaken further the impact of dissipation in 
the metal. 

These remarks are relevant to works suggesting 
the possibility to stabilize the superconductivity in 
ID systems by cou pling them to a bath of normal 
quasiparticle4^^'24125]^ jj^ these works, the basic under- 
lying physical idea is that the normal bath provides a 
source of friction for the phase field 6 (x) which tends 
to quench its fluctuations and therefore, to favor super- 



12 



conductivity (very much like in the case of a resistively 
shunted Josephson junction^ ^^). However, little atten- 
tion has been given up to now to the simultaneous dis- 
sipative effects induced by the Coulomb interaction with 
the electrons in the metal, which produce friction in the 
dual field ps (x), and therefore tends to increase phase 
fluctuations, thus deteriorating superconducting proper- 
ties. In that sense, our results show that the best condi- 
tion would be to screen the Coulomb interaction with 
a clean (i.e., large diffusion constant D) metallic fllm 
(rather than a wire). This result lends credence to the 
analysis made in Ref. [25J where it was assumed that 
the Coulomb interactions only renormalize the bare Lut- 
tinger parameters of a superconducting wire in contact 
with a 2D normal diffusive metal system. 

Many other issues remain to be addressed to get an 
accurate physical description of a superconducting wire 
coupled to a dissipative electron gas, such as the afore- 
mentioned effect of topological excitations^^ , localization 
effects in the gates as a consequence of disorder, simulta- 
neous effect of Coulomb interactions and Andreev tunnel- 



ing, etc. We expect that our results inspire other works 
along these lines. 



Acknowledgments 

This work was supported by the Swiss National Foun- 
dation under MaNEP and division II. 



Appendix A: Derivation of the effective action 

Although the derivation of the low-energy action 
for a supercond uctor has been studied by several 
authord^^ESIllllH here we follow more closely the deriva- 
tion of De Palo et alW^. Our starting point is the decou- 
pling of the interaction terms appearing in and i7i„t 
[Eqs. ([3| and (|7| respectively] by the means of Hubbard- 
Stratonovich transformations (HSTs) 



- Jdr H^„t{r) 



(Al) 



I? [A*, A] e 



' S d'^x^ '^-^^^f^^+J d^x^ A'(x^)ili_l(x^)ili^{x^)+tlj\(x^)xl)l{x^)A{x^) 



(A2) 



where we have introduced the compact notation = 
(r, r) and the bosonic flelds pi, (x^) , A* (x^) , A (x^). 

The quantity [v^ (x^ — a;^)] ^ is a compact notation for 
the Fourier transform 



[v^{x^-x'^)] 



1 ^» 



(A3) 



where (A;^) = ^ d'^x^e ^''''(^f (x^ - a;^) , 

with Vi, {x^ — xjj) the potential (a;^ — x'^ = 



Vi,[r-v')5{T~T'). Note that the mode fc^ = 0, for 
which the above HST is formally ill-defined, can be safely 
ignored by considering the interaction with the positive 
ionic background in the system (not explicitly written 
here). 

Our next step is to decouple the quadratic term 
Pv {Xfi) Pv (^^i) ill Eq. (All by the means of an extra 
HST. According to Ref. \l\A this has the advantage of 
introducing the physical densities (symmetric and anti- 
symmetric) of the problem (cf. Eq. |6|. Then, 



-h f d'^x^d'^x'^ [v„{x^-x'^)] ^ Pu(x^)pu{x\^) 



-\ J d'^x^d'^x'^ Pu{x^)v^[x^,-x'^)p^[x'^)-i J d'^x^ p^(x^)p^(x^,) 



(A4) 



Note that the formal integration of the field pi, gives the 
functional-delta function"^" 5 [p^ — p^] . This fact allows 
to interpret the HS fields pi, as the physical electronic 
densitiea^Sl. 



It is convenient to write the action of the system after 
these manipulations 
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+ y" ^ I??; (9, - M + ^i) + ^ [Vr?:] [V77,] I ^ 

+ j d'^Xfj, pv{x^)[pi,{x^t) - Pv{xit)], (A5) 

where for simplicity we have dropped the arguments in 
the fermionic fields V'o- and rj^ and in the disorder poten- 
tial Vi = Vi{r). 

The next step is to perform the saddle-point 
approximation with respect to the bosonic fields 
A {x^) , A* {x^) , (x^) , pi, (x^), which gives the equa- 
tions 



SS 



_ ^ (^m) 



6A* ix,) \U\ 

= = 



(A6) 



5A (a;^) 

Spi, {x^) 

SS 
Spi, ixf_,) 



\U\ 

= = p^ {Xf_,) - Py (Xf,) , 



(A8) 



= I dx'^ Vy {x^ - x'^) p^ {x'^) + ip^ (x^) . 

(A9) 



The first two equations reproduce the well-known BCS 
gap-equatioEp2lj while the other two give the relation- 
ship between p,^,Pi, and the electronic density. These 
equations provide the starting point for a controlled ex- 
pansion in terms of Gaussian fluctuations of the bosonic 
fields around the uniform solutions A*^"-*, p^^ and pl^\ In 
what follows, we assume that the values of A^°\ and 



are known. When these solutions are inserted back 



into the action Eq. ( A5 I , we notice that the quantity 



pi^'' = ipl?*' J d'^x^ Vy (Xfi) can be absorbed in a renor- 
malization of the chemical potential p^, due to the effect 
of Coulomb interactions, while the divergent quantity 

i J2u=± ('°^"'') / d'^'Xf^d^x'^ Vy [x^ - x'^) exactly cancels 
the contribution coming from the positive ionic back- 
ground (which we have not written explicitly here), by 
imposing the overall electroneutrality of the system, and 
consequently we will drop it in the following. We also 

drop the constant term /Sfiyj^, where Vt is the volume of 
the superconducting system. 

At sufficiently low energies, amplitude fluctuations of 
the order parameter can be neglected, and we can write 
A(a;^) = Aoe'^^^^"), with a real constant Aq = |A(°)|. 



We can absorbe the phase field by the means of a trans- 
formation of the fermion field 

The expression of the effective action is considerably 
simplified introducing the Nambu notation 



■nixf,) 



which allows to write the action as 



Sc^J d^x^ [Ao,s - S,,] ^ + r]^ [Ao,n - S„] v} + 
+ \^ j d^x^d'^x'^ Spy {x^) Vy {x^, - x'^) Spy (x^) + 



+ d*X^Spy{Xf,) p'-y'^ + Spy {Xf,) 



where 



Spy (Xf,) = py (X^) - Pl 



Spy (Xf,) = Py (Xf,) - pI 



(0) 



(AlO) 
(All) 



are the fluctuations of the density around the saddle- 
point solutions, and 



2m 



Ps > T3 - AqTi, 



i{dr9) [vey U «(V6i)v. 



8m 



2m 



where r,; are the Pauli matrices and where we have used 

the fact that ^5 ~ in a translationally in- 

variant system. 

The next step consists in using the expansion formula 

00 / ixn 

Tr In [A) - S] = Tr In A) - V ^Tr [GqS]" , 

^ — ' n 

?i=i 

where Gq = — [.4fl] ^ . Truncating the series at second 
order (i.e., Gaussian fluctuations), we obtain 

S ~ -Tr [Go,.S,] + ^Tr [Go,.£.]' + 

- Tr [Go,„S„] + ^Tr [Go,„S„]' + 

+ ^ X! / d'^Xf.d'^x'^ Spy (Xf,) Vy {x^ - x'^) Spy (x^) + 



+ / d*X^ Spy {Xf,) p^y^+Spy{Xf_,) 

u=±'' 
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where the propagators in Nambu space Go,s and Go.; 
write 



9o,s {x,,) fo,s {x^) 

/0,s (x^) 50,s (Xf,) 



Q ^ / 90,n {Xf,) 

V So.nCa;^) 



The final step is to integrate out the modes dpi, (a;^). 
To that aim, we decouple the mixed term ~ {Spi,) {Spv)' 
appearing in Tr [Go, sSs] and Tr [Go,nSn] by returning 
to the original representation for the densities [cf. Eq. 

(lell 



and where go.sixf,) = - {Tr^^ (Xf,) t/j^ {0)j and 

go,siXf,) = ^T^V'I (a;^) V"; (0)^ denote respectively the 
particle and hole propagators in the superconduc- 
tor, while fo,s{Xf,) = {Tril;i{x^)tp^ {0)), fo,s{Xf,) = 

(^Tt-V'I (xp) i/'l (0)^ are the anomalous ones*'^. Simi- 
larly 5o,n(a;^) = - (TrT]^ {x^)ri^ (0)'^ and .go,„ (a;^) = 

(Tr'T]'^ (x^) 7]^ (0)^ are the particle and hole propagators 
in the normal metal, respectively. 
The evaluation of the traces yields 



Spn 



6p+ + 6p- 
2 

Sp+ - 6p- 



and integrate out the fields Spn (x^) and Sps (x^) instead. 
From here we see that the term ~ {dr9) (drO)' cancels, 
as in Ref. |20l We finally obtain 



S = ^ I d Xf^ drO {x^) ps {x^) 



1 



+ I d^x^d'x'^ \ -V [x^ - x'^) ve (x^) ve {x'^) + 



Tr [Go,.S,] = j d^x^ 

Tr [Go^sS^]^ = j d'^Xf.d'^x'^ |xo,s {x^, - x'^) x 



i^5p^ 



2 ' ^ 2 

^ J Xfj, 



^6p^ {x^,)V {x^,~x'^6p{x'^^ 
^ where we have defined 



(AM) 



(A15) 



Spn (a;^) 



Tr [Go,„S„] = zpio) / d^x,, ^ {i,) Sp, 

If 

Tr [Go,„S„]^ = J d'^Xf_i,d*x'^ xo,n {xf^ ~ x'^) 



Vixy) 



_ / [XO,s {Xy)Y 







■Sir) 



[XO,n (Xy)] 

v{xy,0) v{x^,d) \ 
V {Xfj,,d) V (x^,0) J ' 



^P'^ 



Y,{v')5~p,, 



where again we used a compact notation of Eq. (A3 1. 



where for simplicity we have used the campact notation 
[■^x^ = {xti)^ and where we have defined 

V (x^ - x'^) = [Vgo,s {x'^ - x^) V.go,s (a;^ - a;^,) + 

+V/o.s (a;^ - a;^) V/o,s (a;^ - x^J] , 
and the density-density correlation functions 

XO,s (Xy - X^) = - (TrSps (Xf,) Sps (x^)) , 

= -2go,s (a:^ - x^) 5o,s (a;^ - x^) + 
+ "2/0,8 {xf, - x^) /o,s (x^ - x^) , (A12) 

XO,n {Xy - X^,) = - (TrSp,, (x^) Sp,, (x^)) , 

= -2.go,„ (a;^ - a;^) go,n (a;^ - a:,.) • 

(A13) 



Appendix B: Density susceptibility and 
superconducting stiffness in the limit q ^ 



From Eq. (A12l, the Fourier transforms reads 

xo,. (qn = E [-50,. {kn 9o,s [k^ - qn + 

+/o,s (fc^) /o,. [k-q^]- 
In the limit — > 0, we obtain 



lim xo,. (g^) 



2^ 1 ^ KO'+?k- 



" ^{ii^nf - E'] 



2 

n 



E 



2E. 



2Eu 
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with a = Is - Ais, and E], = + K At T = 



where 7 = f^" = 



2 In 



2 



^2i?w 



(Bl) 



= In 



and where ud is a high-energy cutoff. 

Similarly, the Fourier transform of the superconducting 
stiffness reads 



Ao 



Evaluating the Matsubara sum over the fermionic fre- 
quencies ivm we obtain the result in the limit ^ 



lim V (qn 



-Ek-q 



V ^ 2m2 Ex^ 

k 



1 dnpiEk) 

k 



(0) 
4to 



with 



X [f{k^)f{k^-q^')+g{k^')g{k^'-q' 



-k. (k-q) 

2m2 



^ ]_ ^ A§ + (Wn + gk) (^i^K - ^^m + Ck-q) 



{iVn - iuJmY 



which vanishes in th e limit T 0, and we recover the 
well-know result!221£2l 



JO) 



lim V{qn^Vo = ^. 
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